Ground state and elementary excitations of single and binary 
Bose-Einstein condensates of trapped dipolar gases 
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We analyze the ground-state properties and the excitation spectrum of Bose-Einstein condensates 
of trapped dipolar particles. First, we consider the case of a single-component polarized dipolar gas. 
For this case we discuss the influence of the trapping geometry on the stability of the condensate as 
well as the effects of the dipole-dipole interaction on the excitation spectrum. We discuss also the 
ground state and excitations of a gas composed of two antiparallel dipolar components. 



I. INTRODUCTION 

The nature and stability of a Bose-Einstein condensate 
(BEC) 0J|] are strongly influenced by the interparticle 
interactions, which are described by the s-wave scatter- 
ing length a. If a > the interactions are repulsive, 
and condensates with an arbitrary large number of par- 
ticles are stable. On the contrary, spatially homogeneous 
condensates with a < are absolutely unstable with re- 
gard to local collapses ||. The presence of a trapping 
potential changes the situation drastically, as revealed in 
successful experiments with magnetically trapped atomic 
7 Li (a = —14 A) As found in theoretical studies 

there will be a metastable BEC if the number of 
condensed particles is sufficiently small, such that the 
spacing between the trap levels exceeds the mean field 
interparticle interaction no\g\ (where no is the conden- 
sate density, g — Anft^a/m, and m is the atom mass). In 
other words, the BEC is stabilized if the negative pressure 
caused by the interparticle attraction is compensated by 
the quantum pressure imposed by the trapping potential. 

The effects of the interparticle interactions on the con- 
densate properties have been mainly discussed for the 
case of van der Waals (short-range) interactions. How- 
ever, the BEC in the presence of dipole-dipole interac- 
tions has recently raised a considerable interest [p|-^6[. 
Novel physics is expected for dipolar BEC, since the 
dipole-dipole interactions are long-range, anisotropic and 
partially attractive. The non-trivial task of achieving and 
controlling dipolar BECs is thus particularly challenging. 

The interest on dipolar gases has been partially moti- 
vated by the recent success in creating ultra-cold molec- 
ular clouds p7 18 1 . This success opens fascinating 
prospects to achieve quantum degeneracy in trapped 
gases of heteronuclear molecules, which could interact 
via electric dipole-dipole forces after being oriented in a 
sufficiently high electric field. On the other hand, the 
ultra-cold gases of atoms with large magnetic dipole mo- 
ments, as chromium [[l9|-p3[ or europium p4[ , have also 
been subject of growing interest. In this case, the dipole- 
dipole interactions are not expected to be dominant, al- 



though for a relatively small a the BEC may reflect the 
interplay between short-range and dipole-dipole interac- 
tions pjlgf|. Interestingly, these effects can be amplified 
by reducing a via Feshbach resonanc es p5l p6||. Similar 
effects have been discussed in Refs. |g|, |ll| , |12f for ground- 
state atoms with electric dipole moments induced by a 
high dc electric fields (of the order of 10 6 V/cm). It 
has been also suggested that laser-induced dipole-dipole 
interactions could be achieved by exciting atoms to Ry- 
dberg states pX[ |. In this case applications to quantum 
information processing have been discussed (see e.g. p7j). 

The stability of the condensate is significantly modified 
by the presence of dipole-dipole interactions [p[- pl]JT^ ] . In 
particular, a BEC of particles dominantly interacting via 
dipole forces is, similarly to condensates with a < 0, un- 
stable in a spatially homogeneous case and can be sta- 
bilized by confinement in a trap. It has been shown |lC| ] 
that the sign and the value of the dipole-dipole inter- 
action energy is strongly influenced by the trapping ge- 
ometry and, hence, the stability diagram depends cru- 
cially on the trap anisotropy. This opens new interesting 
possibilities for controlling and engineering macroscopic 
quantum states. In particular, for dipoles oriented along 
the axis of a cylindrical trap there exists a critical value 
U = 0.4 for the square root of the ratio of the radial to ax- 
ial frequency I = (ojp/ujz) 1 / 2 . Pancake traps with I < I* 
provide mostly a repulsive mean field of the dipole-dipole 
interaction, and thus the dipolar condensate in such traps 
will be stable at any number of particles N. For / > Z* 
the stability requires N < N c , where the critical value N c 
at which the collapse occurs is determined by the condi- 
tion that (on average) the mean field interaction energy 
is attractive and close in absolute value to hui p . 

The study of the condensate properties would be in- 
complete without the analysis of the excitation spectrum, 
which determines the dynamical behavior of the system 
in the regime of weak perturbations. Zero-temperature 
excitation frequencies have been extensively studied in 
the case of condensates in dilute alkali gases, both exper- 
imentally |2§| , |29f| and theoretically f30|^32||. In this paper 
we discuss the low-lying collective excitation frequencies 
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of trapped dipolar BEC, complementing the analysis of 
Refs. OJI3| . We first consider the case of dominant 
dipole-dipole interactions, and later on we discuss the 
situation where the short-range interaction is also rele- 
vant. In particular, we discuss in detail the nature of 
instability and demonstrate that one of the excitations 
frequencies tends to zero at the criticality as a power of 
(N c — TV)' 3 . We discuss the qualitative character of the 
low- lying modes, and show that the exponent (3 under- 
goes a crossover from 1/4 for I ^> Z» to 2 at I ~ Z». To 
illustrate the case of mixed short-range and dipole-dipole 
interactions, we present predictions for the excitation fre- 
quencies for the particular case of a chromium gas. 

The second part of the paper is devoted to the analysis 
of binary dipolar BECs. In recent years the development 
of trapping techniques has allowed for creation of mul- 
ticomponent condensates, formed by atoms in different 
internal (electronic) states |j3|]34| ]. The physics of mul- 
ticomponent BEC, far from being a trivial extension of 
the single-component one, presents novel and fundamen- 
tally different scenarios for its ground-state wave func- 
tion |35],[36| and excitations |37]]. In particular, it has 
been experimentally observed that a BEC can reach an 
equilibrium state characterized by phase separation of 
the species in different domains Hj. The analysis of 
multicomponent BEC has been so far mostly limited to 
the case of short-range interparticle interactions (a model 
long-range interaction has been considered in Ref. [p6| ). 
One of the main goals of this paper is to analyze the prop- 
erties of multicomponent dipolar BEC. Such a mixture 
can be achieved in different physical contexts. In par- 
ticular, it would be in general the case for experiments 
in ultracold polar molecules |58|, and in chromium ]39| ] 
in which different magnetic-moment species are simulta- 
neously trapped in an optical dipole trap. It would also 
be the case of atomic electric dipoles created by laser- 
induced pumping to two different Rydberg states. Fi- 
nally, the same situation would appear in condensates of 
heteronuclear Hund A diatomic molecules, for which the 
direction of the magnetic moment is correlated (parallel 
or antiparallel) with the direction of the molecular axis. 
Thus, if the magnetic moments are oriented in a mag- 
netic field, the electric moments can acquire two possible 
directions. We show below that a binary dipolar BEC of 
two antiparallel dipole components differs qualitatively 
from the case of a short-range interacting binary BEC. 

Our paper is organized as follows. In Sec. || we briefly 
review the ground-state properties of a single component 
BEC of trapped dipolar gases JTc|] . Sec. |n| is devoted 
to the analysis of the excitation spectrum of single com- 
ponent trapped dipolar condensed gases. Sec. IV briefly 
discusses the ballistic expansion of a dipolar BEC. In 
Sec. [y| the ground state of a BEC of two different dipolar 
species is considered. The excitation spectrum for this 
case is discussed in Sec. VI. We conclude in Sec. VII. 



II. GROUND-STATE PROPERTIES OF A 
SINGLE-COMPONENT DIPOLAR BEC 

A. Description of the system 

In this section we briefly review the results of Ref. 
|l0| . We consider a condensate of dipolar particles in 
a cylindrical harmonic trap. In the following we consider 
the case of electric dipoles, although the same physics 
is expected for magnetic ones p. All dipoles are as- 
sumed to be oriented along the trap axis by a suffi- 
ciently strong external field. Accordingly, the dipole- 
dipole interaction potential between two dipoles is given 
by Vd(R) = (d 2 /R 3 )(l — 3cos 2 #), where d characterizes 
the dipole moment, R the vector between the dipoles 
(R = |R| being its length), and 9 the angle between R 
and the dipole orientation. Similarly as in Refs. | 0, 
we describe the dynamics of the condensate wave func- 
tion ^(r, t) by using the time-dependent Gross-Pitaevskii 
equation (CPE): 
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hKr',t)| 2 U(r,i). (1) 



where ip(T,t) is normalized to the total number of con- 
densate particles N. The third term in the rhs of Eq. 
(|l|) corresponds to the mean field of short-range (van 
der Waals) forces and the last term to the mean field 
of the dipole-dipole interaction. In this section, we 
omit the term g\ifj(r,t)\ 2 ip(r,t), assuming that the inter- 
particle interaction is dominated by dipole-dipole forces 
(d 2 3> \g\ = AnTi 2 \a\/m) and the system is away from 
shape resonances of Vd (R) . The effects of the short-range 
interactions on the excitation spectrum are discussed in 



detail in Sec. Ill 



The wave function of the relative motion of a pair of 
dipoles is influenced by the dipole-dipole interaction at 
interparticle distances |r — r'| < r* = 2md 2 /h 2 . This 
influence is ignored in the dipole-dipole term of Eq. ([!]), 
as the main contribution to the integral comes from dis- 
tances |r — r'| of order the spatial size of the condensate, 
which we assume to be much larger than r*. 

Away from the shape resonances the dipolar conden- 
sate is unstable in the spatially homogeneous case. For 
all dipoles parallel to each other, by using the Bogoliubov 
method, one finds an anisotropic dispersion law for ele- 
mentary excitations: e(k) = [E 2 + %tt Ekn^d 2 {cos 2 9 k — 
1/3)] 1 / 2 , where E k = h 2 k 2 /2m, no is the condensate 
density, and 9^ is the angle between the excitation mo- 
mentum k and the direction of the dipoles. The insta- 
bility is clearly seen from the fact that at small k and 
cos 2 9k < 1/3 imaginary excitation energies e emerge. 
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B. Numerical results 

Equation (|IJ), contrary to the usually employed GPE 
with short-range interactions, is an integro-differential 
equation. The evaluation of the integral term deserves 
special attention, since the integrand diverges at rela- 
tive interparticle distances tending to zero. Fortunately, 
the calculation of the integral term can be simplified by 
means of the convolution theorem: 



dr 



,1-3 cos 2 

- r-'|3 



#(rO| 2 = ^-M^](q»l 2 ](q)}, 

(2) 



where T and T~ x indicate the Fourier transform and 
the inverse Fourier transform, respectively. The Fourier 
transform of the dipole-dipole potential reads: 



F[V]{q) = 4nd 2 {l - 3 cos 2 a) [ 



cos(bq) sin(bq) 



(bqf 



(bqf 



(3) 



where a is the angle between the momentum q and the 
dipole direction, and b is a cutoff distance corresponding 
to the atomic radius (a few Bohr radii). Since b is much 
smaller than any significant length scale of the system, 
one can safely perform the limit 



lim T(V(r)) 



in 



-d 2 (3 cos 2 



1). 



(4) 



In order to evaluate the Fourier transform of Eq. (g) 
J r (|V ; | 2 ) is numerically evaluated by means of a standard 
Fast Fourier Transform (FFT) algorithm and multiplied 
by T(V(y)). 

The ground-state of the system is obtained by em- 
ploying a standard split-operator technique in imaginary 
time. The split-operator technique is also based in an 
FFT algorithm and, consequently, for each time step four 
FFT's are needed: two for the calculation of the integral 
term and two for the evolution. We would like to stress 
that this procedure constitutes a non trivial computa- 
tional task. Additionally, the FFT algorithm must be 
evaluated in Cartesian coordinates and, as a consequence, 
computationally demanding fully three-dimensional cal- 
culations are required. 

To understand the influence of the trapping potential 
on the dipolar condensate, we have simulated Eq.(|l|) for 
various values of the number of particles N, dipole mo- 
ment d, and the trap aspect ratio I. We have found the 
conditions under which the condensate is stabilized by 
the trapping field and investigated static properties of 
this Bose-condensed state. 

For a stationary condensate the wave function 
ip(r, t) = ipo(r) exp (—ifit/h), where /i is the chemi- 
cal potential, and the lhs of Eq. ([j]) becomes [iipo(r). 
The important energy scales of the problem are the 
trap frequencies uj z , oj p and the dipole-dipole interac- 
tion energy per particle defined as V — (1/N) J Vd(r — 



r')ip 2 (v)ip^(r')drdr' . Accordingly, the trap frequen- 
cies, and the (renormalized) number of particles a = 
Ar»/a max (with a max = (h/2mui m i n ) 1 ^ 2 being the maxi- 
mal oscillator length of the trap) form the necessary set 
of parameters allowing us to determine the chemical po- 
tential and give a full description of the ground state of 
a trapped dipolar condensate. 

We have found that the dipolar condensate is stable 
either at V > or at V < with |V| < Tiu p . This re- 
quires N < N c , where the critical number N c depends on 
the trap aspect ratio I = (ujp/uJz) 1 ^ 2 . The calculated de- 
pendence N c (l) clearly indicates the presence of a critical 



point Z* ~ 0.43 [ 
condensate is sta 
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In pancake traps with I < Z* the 
Die at any TV, because V always remains 
positive. For small TV the shape of the cloud is Gaus- 
sian in all directions. With increasing TV, the quantity 
V increases and the cloud first becomes Thomas-Fermi 
in the radial direction and then, for a very large TV, also 
axially. The ratio of the axial to radial size of the cloud, 
L = L z /L p , continuously decreases with increasing num- 
ber of particles and reaches a limiting value at TV — > oo 
(see Fig. 3 of Ref. [HJ). In this respect, for a very large 
TV we have a pancake Thomas-Fermi condensate. 

For I > 1 the mean field dipole-dipole interaction is al- 
ways attractive. The quantity |V| increases with TV" and 
the shape of the cloud changes. In spherical traps the 
cloud becomes more elongated in the axial direction and 
near N = N c the shape of the cloud is close to Gaus- 
sian with the aspect ratio L — 2.1. In cigar-shaped 
traps (! > 1) especially interesting is the regime where 
huj z <C |V| <C fiujp. In this case the radial shape of the 
cloud remains the same Gaussian as in a non-interacting 
gas, but the axial behavior of the condensate will be 
governed by the dipole-dipole interaction which acquires 
a quasi 1- dimensional (ID) character. Thus, one has a 
(quasi) ID gas with attractive interparticle interactions 
and is dealing with a stable (bright) soliton-like conden- 
sate where attractive forces are compensated by the ki- 
netic energy po| . With increasing TV, L z decreases. Near 
TV = TV C , where \V\ is close to hw p , the axial shape of the 
cloud also becomes Gaussian and the aspect ratio takes 
the value L w 3.0. For Z* < I < 1 the dipole-dipole inter- 
action energy is positive for a small number of particles 
and increases with N . The quantity V reaches its maxi- 
mum and further increase in N reduces V and makes the 
cloud less pancake-shaped. At the critical point N = N c 
the shape of the cloud is close to Gaussian and the aspect 
ratio L < 3.0. 

In the previous analysis, the case of dominant dipole- 
dipole interactions was considered. However , 
in the general case the effects of the short-range inter- 
actions are comparable or even larger than those related 
to the dipole-dipole interactions. In such situations the 
short-range term must be maintained in Eq. (]]]). Pro- 
vided that a is sufficiently small and positive, the system 
can become unstable and undergo a collapse in a similar 
way to what was observed in experiments with 7 Li 
and 85 Rb Efl]. For the case of negative a the dipolar gas 
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is expected to be highly unstable, but the dipole-dipole 
interaction could be employed to stabilize the gas by the 
trap geometry in a way analogous to the one presented 
in Ref. fujl. 



III. EXCITATIONS OF A SINGLE-COMPONENT 
DIPOLAR BEC 

In this section we analyze the collective excitations of a 
dipolar BEC. Since this is a potentially unstable system, 
there is a fundamental question about the qualitative and 
quantitative nature of this instability. Another impor- 
tant question concerns quantitative aspects, in particular 
how relevant are the effects of the dipole-dipole interac- 
tion in the excitation spectrum and to what extent they 
can observed in the exper iments. 

In the subsection III A , we briefly summarize the vari- 
ational approach introduced by Perez-Garcfa et al. ]3l| ] 
and later used by Yi and You |ll]Jl2]| to describe the low- 
lying excitation spectrum of a dipolar condensate (for a 
different variational method using the self-similarity as- 
sumption see |41|). We employ t he no tation introduced 
in the Refs. |11| , |12| . Subsection III B| is devoted to the 
analysis of the behavior of th e excitations close to the 
instability. In subsection IIIC, we calculate numerically 



the response of the system to small perturbations and 
compare the result s with those yielded b y the method 
of subsection 



IIIB 



Finally, in subsection III C , we also 
discuss how the effects of the dipole-dipole interactions 
can manifest themselves in the excitation spectrum for 
various kinds of dipolar gases and trapping geometries. 



A. Variational approach 

The problem of solving the time dependent GPE Eq. 
(|l|) can be restated as a variational problem correspond- 
ing to the stationary point of the action S = j dtdvC 
related to the Lagrangian density 0,n3,Rl[ : 



£ = -h 
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r — r 



(5) 



where u-i = i^/nJEd 2 and Vt(r) describes the trapping 
potential. We consider the following Gaussian ansatz for 
the condensate wave function: 



1>(x,y,z,t)=A(t) J] e -'7 2 /2^= e in 2 ^(t) ; 



(G) 



ansatz (||) into Eq. (||) and integrating over the spa- 
tial coordinates, one obtains the effective Lagrangian. 
From this effective Lagrangian one finds the equations 
of motion of the variational parameters, i.e. the cor- 
responding Eulcr-Lagrange equations. In particular, if 
w n = v^^h/muj, LO = U0 x . y , L0 Z = \u x , v XiV = V, 

p 



T 



/2/irNa/ \Jh/mu}, and r = tot, the Euler-Lagrange 
equations take the following form: 
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drexp — 



V- If 



Y w (0) 



where X x ^ y = 1, X z = A. This equation describes the 
motion of a particle with coordinates (v x , v y , v z ) in an 
effective potential 
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- K4 + ^z v z) 
1 1 \ p 



V X VyV Z 



dr exp 




mode 1 



mode 2 



mode 3 



(7) 



(8) 



FIG. 1. Graphical representation of oscillations modes of 
the condensate. 

Small amplitude oscillations around the stationary so- 
lution can be found by studying second derivatives of U . 
This procedure provides the excitation frequency for the 
first three compressional excitation modes of the system 
pl| . The geometry of these modes is depicted in Fig. ||. 
For cylindrical traps with the axis along the dipole direc- 
tion, the projection of the angular momentum, to, on the 
z axis is a good quantum number. For modes 2 and 3 we 
have to = 0, whereas to = 1 for mode I. In the following 
we call mode 2 (3) as the breathing (quadrupole) mode. 



B. The nature of the instability 



where in contrast to Ref. [jlTJ we omit the possibility of 
the sloshing motion of the condensate. By inserting the 



To study the nature of the instability one needs to de- 
termine which mode becomes unstable when the number 
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of atoms reaches the critical value and to describe the be- 
havior of the frequency of the corresponding mode close 
to the criticality. These two issues have been first dis- 
cussed for the case of short-range interacting Bose gases 
with a < 0. Bergeman (32] observed numerically that, as 
the ratio 7 of the nonlinear interaction energy to the trap 
frequency approaches a given critical value 7 C , the fre- 
quency of the breathing mode tends to zero and merges 
with the frequency of the Goldstone mode correspond- 
ing to the overall phase of the condensate. Above the 
criticality, the breathing mode becomes unstable and at- 
tains complex frequency. Singh and Rokhsar ]30| | ana- 
lyzed this instability using self-similar solutions describ- 
ing the modes (equivalently one may employ the varia- 
tional approach of the previous subsection). They have 
shown that close to the criticality the frequency of the 
breathing mode 2 vanishes as I7 — 7 C | 1 ^ 4 - 

In the case of a dipolar gas with dominant dipole in- 
teractions the situation is completely different and much 
more complex. Only for aspect ratios far above the crit- 
icality, I ^S> I* (I > 1.29) the situation resembles that of 
a gas with a < 0. The mode corresponding to the lowest 
frequency is the breathing mode. This mode becomes un- 
stable when the parameter a — > a c . The scaling behavior 
of the frequency of this mode can be analyzed employing 
the variational approach of the previous section 11 1 . We 



V(r,t) = ^m 

7j=x,y,z 



[1 + A v sin(u mod t + a v )] 2 oj%T) 2 , (9) 



find that u>2 goes to zero as (<r c — a)", with f3 ~ 1/4. 
For intermediate values of I > I* (0.75 < I < 1.29) the 
exponent (3 is still 1/4, but the geometry of the mode 
which drives the instability depends on a. For a far be- 
low cr c the mode corresponding to the lowest frequency 
has a breathing symmetry, whereas as one approaches the 
critical value of a the lowest mode becomes quadrupole- 
like. For I close to I* (I < 0.75) the situation changes 
and the mode corresponding to the lowest frequency is 
quadrupolar-like. For I not too close to 1* the exponent /3 
is still close to 1/4. However, as / approaches I*, the ex- 
ponent [3 departs from 1 /4 towards a greater value, which 
becomes 2 at I = ?*. The latter prediction must be un- 
derstood as a qualitative argument, since for I very close 
to the system enters a Thomas- Fermi regime where 
the Gaussian ansatz is no more valid. 



C. Numerical results 

In this subsection we study the low-lying collective ex- 
citations of a trapped dipolar BEC by analyzing numer- 
ically the response of a BEC after applying an external 
perturbation. We compare our numerical results with 
those obtained from the variational approach described 
above. In the following, we employ the method intro- 
duced in Ref. 0] for the case of a short-range interact- 
ing BEC. After generating the BEC ground-state wave 
function using imaginary time evolution of the GPE, we 
first slightly perturb the trapping potential in a periodic 
way: 



where Lu mo d is the modulation frequency, A v are the am- 
plitudes, and otrj are the initial phases. A perturbation of 
a chosen symmetry can be accomplished by properly se- 
lecting these parameters. In particular, sufficiently small 
amplitudes are necessary near the instabilities, since the 
system can be easily driven into collapse. On the other 
hand, by setting large amplitudes one may probe various 
nonlinear effects [^2|,£l3| . The response of the system is 
enhanced by selecting oj mo d in the vicinity of expected 
excitation frequencies. 

In a second stage, the condensate evolves in an unper- 
turbed trap (i.e. A v — 0). The condensate widths are 
monitored and subsequently Fourier-transformed to re- 
veal the excitation spectrum. Our aim is to determine 
the (potentially small) deviations of the excitation fre- 
quencies with respect to those expected for purely short- 
range interacting gases, as precisely as they can be exper- 



imentally measured (typically Auj/cotrap — 0.01 



)■ 



Long integration times are necessary to accomplish the 
desired spectral resolution. This introduces serious tech- 
nical difficulties, since a large number of integration steps 
is needed to guarantee the energy conservation during the 
whole evolution. Note that additional complications arise 
from the evaluation of the interaction integral in Eq. (Q) 
at each time step. 



1. Dominant dipole- dipole interactions 

For the case of dominant dipole-dipole interactions 
JTo| | , the short-range part of Eq. ([!]) can be safely omit- 
ted. We consider two different trap geometries, that were 
employed at JILA in Ref. p8| , with a trap aspect ratio 
I = (referred to as pancake JILA trap), and a trap 

with I — 8 1 / 4 (referred to as cigar JILA trap). The later 
trap has been recently employed for experiments in 85 Rb 

m 

In the following, we study the excitation frequencies as 
a function of the dipolar parameter 



( = N- 



(10) 



Note the relation a — 2y / '2ujp/uj m i n (. Since in the pan- 
cake JILA trap the breathing mode has a relatively high 
frequency and it is not excited by the perturbation of 
the GPE, we concentrate in this case on the modes 1 
and 3 of Fig. [l[ The results for the pancake JILA trap 
are presented in Fig. |[ For £ = the ideal-gas result 
lu/lo p — 2 is retrieved for both modes 1 and 3. As £ in- 
creases mode 1 is shifted upwards, whereas the frequency 
of mode 3 (quadrupole) goes down, vanishing for a crit- 
ical £ cr at which the system becomes unstable. Probing 
the system in this regime is very difficult, since even a 
slight disturbance drives the collapse. For £ > ( cr ss 3.68 
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it is not possible to obtain stable ground state solutions 
of the time-independent GPE |l(J . The variational anal- 
ysis reproduces the numerical results for relatively weak 
dipolar interactions. However, it does not describe well 
the numerical results close to the instability. 
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FIG. 2. Numerical results for the excitation frequencies of 
modes 1 (squares) and 3 (circles) as functions of the dipolar 
parameter £ for a gas of condensed dipoles in the pancake 
JILA trap. The solid lines indicate the corresponding varia- 
tional results. 

Fig. H depicts a similar dependence for the cigar JILA 
trap. In this case the breathing mode is the lowest 
one. The frequency of mode 1 displays an upward shift 
whereas that of the quadrupole mode stays essentially 
untouched. Again, the instability threshold is not cor- 
rectly predicted by the variational analysis. 




0.5- 



FIG. 3. Numerical results for the excitation frequencies of 
modes 1 (squares) , 2 (triangles) and 3 (circles) as functions of 
the dipolar parameter £ for a gas of condensed dipoles in the 
cigar JILA trap. The solid lines indicate the corresponding 
variational results. 



2. Interplay between short-range and dipole-dipole 
interactions 

Let us now consider the case in which the short-range 
interactions and the dipole-dipole ones have a compara- 
ble strength §,|,|l|,^|,|lj . In particular, we have per- 
formed our simulations for the particular case of 52 Cr, 
which has drawn some experimental interest [^| E3|, 
since it possesses a large magnetic dipole moment of 
6/is (Bohr magnetons). The value of a for this ele- 
ment is at the moment unknown, and consequently we 
have explored different values of a in the calculations be- 
low. In the following, the trap parameters correspond to 
those of an ongoing experiment in Stuttgart |39f| , namely 
uj z = 2tt x 40 Hz and uj p = 2tt x 485 Hz. 




FIG. 4. Excitation frequencies of modes 1 (filled squares), 
highest (filled triangles) and lowest (filled circles) as functions 
of the scattering length (expressed in units of the correspond- 
ing value for sodium) for a condensate of 10000 chromium 
atoms in a trap with ui z — 2-7r40 Hz and ui p = 2-7T485 Hz - 
numerical results. The corresponding variational values are 
depicted by dashed lines. Solid lines indicate variational data 
for the case with no dipole-dipole interactions and the corre- 
sponding numerical results are plotted with empty symbols 
(note: we have not been able to probe the lowest mode nu- 
merically) . 

Fig. [| shows the three lowest modes for different val- 
ues of a, for the case of N = 10000 atoms. When a 
is smaller than a critical a cr u the system becomes un- 
stable due to an unbalanced attractive component of 
dipole-dipole interactions [§j^ [13|Jlj] ] . When approach- 
g ing a c "*/°Na ~ 0.23 from above, the frequency of the 
lowest mode decreases to and the system collapses. 
For 0.23 < a/a,Na < 0.26a,/v a the lowest mode (tri- 
angles) has a breathing geometry, whereas the highest 
one (circles) is quadrupole-like. The opposite is true for 
a/iNa > 0.26ciNa- The highest mode shows virtually 
no dependence on the scattering length whereas the fre- 
quency of mode 1 is shifted upwards. Again, close to 
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the instability, the results obtained from the variational 
calculation differ from those obtained numerically. More- 
over, the behavior of the lowest mode is not reproduced 
correctly even qualitatively. Fig. || also presents the exci- 
tation frequencies for the same system in the absence of 
dipole-dipole interactions. Apart from the obvious pres- 
ence of the instability in the dipolar case (and the corre- 
sponding behavior of the lowest mode) we observe large 
frequency differences for mode 1 between the cases with 
and without dipole-dipole interactions. 

The above result suggests that the effects of dipole- 
dipole interactions may be detectable. Having discussed 
the role of a, it is thus important to discuss possible val- 
ues of N and lu p /lu z which could maximize the dipole- 
induccd frequency shifts. 
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FIG. 5. Difference of excitation frequency for mode 1 in the 
absence and in the presence of dipole-dipole interactions, as 
a function of the total nu mber of atoms N, for a — dNa and 
a trap aspect ratio I = U 485/40 (data from the variational 
analysis) . 

Fig. H shows the deviation between the case with 
and without dipole-dipole interactions, for the mode 1, 
lu z /lu p = 40/485, and a = aN a - For this calculation we 
have employed the variational method as it is very exact 
away from the instability. We observe that for number of 
atoms N greater than 5000 the frequency shift is observ- 
able (i.e. greater than 0.01), and reaches its maximum 
(> 3.5%) for N ~ 10000 which should be the number 
available in ongoing experiments p9[| . In this case, shifts 
for modes 2 and 3 are substantially lower (~ 0.001). 

Finally, we have analyzed (using the variational 
method) the dependence of the spectrum on the trap as- 
pect ratio I for a = aNa and N = 10000. The results are 
presented in Fig. ^. We observe that in pancake traps 
with I < 0.7, the quadrupole mode experiences a shift 
of up to 10% with respect to the nondipolar case. We 
have confirmed this conclusion by performing an exact 
numerical simulation of the perturbed GPE for I = 0.58, 
a = fl]\ro, and N — 10000. For these parameters the ex- 
citation frequency of mode 3 is 1.76 to p (1.86 tu p ) with 



(without) dipole-dipole interactions. The remaining two 
modes also display large shifts over a wide range of aspect 
ratios. 
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FIG. 6. Difference of excitation frequency for modes 1 
(dashed line), breathing (dotted line) and quadrupole- like 
(solid line) between the cases of purely contact and mixed 
(contact and dipole-dipole) interactions as a function of the 
trap aspect ratio for 10000 atoms and a = aNa (data from 
variational analysis). 



IV. BALLISTIC EXPANSION OF A DIPOLAR 
BEC 

In typical BEC experiments, the measurements are 
performed after removing the trapping potential and al- 
lowing for a ballistic expansion of the condensate. It is 
the aim of this section to briefly discuss such an expan- 
sion, by means of the previously introduced variational 
approach. After evaluating the minimum of the potential 
provided by Eq. (||), i.e. the ground state of the trapped 
gas, we set \ x , y ,z = in Eqs. (0) and evolve these equa- 
tions using the ground state as an initial condition. In 
this section we restrict ourselves to the case of dominant 
dipole-dipole interactions. 

The ballistic expansion of a short-range interacting 
BEC is characterized by an inversion of the aspect ratio of 
the condensate cloud, i.e. cigar-shaped condensates be- 
come pancake-like after expansion, and vice versa. This 
is not necessarily the case in dipolar condensates. Fig. Q 
shows the condensate aspect ratio for a spherical trap, be- 
fore releasing the trap (dashed), and once the aspect ra- 
tio reaches a stationary value after the expansion (solid). 
The aspect ratio of the cloud decreases during the ex- 
pansion, but never becomes smaller than 1, i.e. the 
condensate keeps its original cigar-shaped character dur- 
ing the expansion. A similar behavior is observed for 
1 < I < 1.12, where there exits a range of £ values for 
which the BEC remains cigar-like. For Z* < I < 1, for 
£ close to ( cr the BEC also remains cigar-shaped during 
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the time of flight. On the contrary, for I > 1.12 the cloud 
becomes pancake-like after expanding for all £, and for 
I < Z» the initial pancake cloud always becomes cigar- 
shaped. 




FIG. 7. Aspect ratio as a function of a for a condensate 
initially trapped in a spherical trap of frequency ui, before the 
expansion (dashed) and after 48.5cj -1 (solid). 



V. GROUND STATE OF A BINARY DIPOLAR 
BEC 

In the previous sections we have shown that various 
novel phenomena are expected in dipolar BEC. One can 
therefore expect that even a richer behavior can be dis- 
played by a multicomponent dipolar BEC. In the case of 
two-component dipolar gases, the system properties de- 
pend on a large number of control parameters, including 
the number of particles in each component, the strengths 
and orientations of the dipoles, and the trap geometries. 
In the following, for simplicity, we focus on the case of 
a BEC of two dipolar identical components polarized in 
opposite directions. As discussed in Sec. | this should be 
the case for Hund A molecules in a magnetic field, where 
the magnetic moment is oriented by the field, but the 
(possibly large) electric moment can be parallel or an- 
tiparallel to the direction of the applied magnetic field. 
To simplify the analysis even further, we consider only 
a spherical trap. Thus, the system is determined by £ 
and rj = Ni/N, where as before N is the total number of 
particles and A^i is the number of particles in the com- 
ponent 1. The system can be described by a system of 
two coupled GPEs: 



V 2 
2 



— T" + -W + +((Wi(r) - (1 - v)V 3 {t)) rh, (Ha) 



2 



+ V + +((-Wi(r) + (1 - v)V 2 (r)) 1P2, (Hb) 



where V^r) = / dr'(l - 3 cos 2 0)|^i(r', i)| 2 /|r - r'| 3 . In 
the above expression we have employed harmonic oscilla- 
tor units of length d = ^h/mw. For equal population of 
the components (77 = 0.5) both ground-state wave func- 
tions are the same. As a consequence, the nonlocal in- 
teraction terms vanish and the system behaves like an 
ideal gas, the ground-state wave functions being that of 
the corresponding harmonic oscillator. However, as the 
interaction strength £ increases, the system becomes un- 
stable. The reason of this instability is rooted in the 
excitation spectrum and is discussed below. 
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FIG. 8. Stability diagram for a binary condensate of oppo- 
sitely polarized dipolar gases in a spherically symmetric trap. 
Above the line there are no stable ground-state solutions. 

For each value of rj it is possible to compute a criti- 
cal £ cr above which no stable solution of Eq. (|ll|) exists. 
The corresponding stability diagram, resulting from a nu- 
merical solution of Eqs. (pT[), is presented in Fig. 0. We 
observe that the more symmetric the mixture is the more 
stable it is. 

For i] 7^ 0.5 the ground state has a more compli- 
cated structure. In this respect it resembles solutions 
for short-range interacting binary condensates p5| , |36|] . 
In particular, the component whose self interaction is 
stronger partially wraps around the other. However, con- 
trary to short-range interacting binary BEC, due to the 
anisotropy of the dipole-dipole interactions, the density 
dip in one component appears only along the radial di- 
rection, whereas along z both components possess Gaus- 
sian profiles. An example of such a structure is shown in 
Fig. I 

In the analysis presented in this section we have ne- 
glected the influence of short-range interactions. Since 
the two components are identical except for the opposite 
dipole moments, the short-range interactions (both inter- 
component and the intracomponent ones) are the same. 
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Therefore, in the absence of dipole-dipole interactions a 
miscible mixture is expected, which differs from the ideal- 
gas Gaussian case considered above (eventually acquiring 
a Thomas- Fermi profile). However, in the same way as 
above, the dipole-dipole mean fields will be exactly can- 
celed for 77 = 0.5. Therefore, the cloud is expected to re- 
main in an unperturbed Thomas- Fermi profile until the £ 
coefficient reaches a critical value £ cr (only quantitatively 
different from that obtained in absence of short-range in- 
teractions). For £ > £ cr the system becomes unstable. 




FIG. 9. Density profiles along the radial (upper graph) and 
axial (lower graph) directions for a binary condensate of two 
oppositely polarized dipolar gases. Parameters are rj — 0.8 
and C = 4.8. 



with 

5(k'i,k2;ki,k 2 ) = 

dr J dr'^W^CrO^Cr-rOVk.CrOV'k^r), (13) 

and aki the annihilation operator of a particle of the i-th 
component with momentum k. Performing the Bogoli- 
ubov approximation (aji = aoi = V^i an d o,q 2 — ^02 = 
\fN~2) : one obtains 

H = ^{[ £ k + 2B(k)]a kl a k i + [e k + 2A 2 S(k)]aj c2 a k2 



-s(k)[(4 



kl"-kl 



akia-ki) + A 2 (ak2 a -k2 + ak2a-k2) 
-2A(a k ia_ k 2 + a ki a-k2 + a k:i a k2 + a k 2 a ki)]}; ( 14 ) 

with ek = fi. 2 k 2 /2m, y^rff = n id\ and -B(k) = 
^ y^(nidi)Y2o(k). The diagonalization of the Hamil- 
tonian provides two different branches of quasiparticlc 
excitations: 



(15a) 




32tt hr 



(l + A 2 )(n 1 d 2 )y 2 o(k)e k . (15b) 



For the case of A — > the spectrum of a homogeneous 
single-component dipolar BEC is discussed in Sec. |J. 
As observed, the spectra (15a) is that of a soft mode, 
whereas the second branch (15b) is even more unstable 
than that of a single component BEC. 



VI. EXCITATIONS OF A BINARY DIPOLAR BEC 



B. Trapped case 



In this section we analyze the excitation spectrum of 
the previously discussed gas of two oppositely oriented 
dipolar components. We first discuss the homogeneous 
case, and then study the case of harmonically trapped 
binary dipolar BEC. 



A. Homogeneous space 

Let us consider two species of dipoles di = d\U z and 
d2 = — d2U z with respective densities n\ and n 2 in a ho- 
mogeneous box of volume V. The Hamiltonian in second 
quantization is: 



2 1,2 



,t 



kct=l,2 

+ \ E E 5 , (ki,k 2 ;ki,k 2 )aj c , Q aJ t , Q a k2Q ak 1 a 

a=l,2ki,k 2 ,ki,k^ 
E 3(ki,k 2 ;ki,k 2 )aJ i , il aj c;)2 ak 2 2ak 1 i, (12) 

ki,k2,k' 1 ,k 2 



For simplicity we limit ourselves to the case of a spher- 
ical trap with frequency lo, but similar arguments can 
be employed for a cylindrical trap. We also consider 
|di| = |d 2 | and 77 = 0.5 (i.e. N x = N 2 = N/2). As we 
have already noted, in this case the ground state of the 
system is the ground state of a harmonic oscillator. It is 
therefore more convenient to rewrite the Hamiltonian in 
terms of the creation and annihilation operators of the 
harmonic oscillator modes for the dipole components 1 
and 2: 

n ni,n 2 #0 

{3(111, n 2 ; 0, Oja^a^ + 5(0, 0; ni, n 2 )a ni ia n2 i 
+fl(ni, n 2 ; 0, 0)a' ni2 a' a32 + g(0, 0; ni, n 2 )a ni2 a n22 
+23(111, 0; n 2 , 0)a| lil a n2 i + 2g(ni, 0; n 2 , 0)a^ i2 o„ 2 2 
-23(0, 0; ni, n 2 )a ni ia n2 2 - 2 5 (ni, n 2 ; 0, Ojo^a^ 
-23(111, 0; 0, n 2 )a] lll a n2 2 - 23(111, 0; 0, n 2 )a| li2 a n2 i}, (16) 
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where n = (n, I, m) denotes the set of spherical quantum 
numbers for the corresponding eigenstate of the harmonic 
oscillator, Kuj n is the energy of the eigenstate n, and 

0(111,112; 113,114) = 

dr J dr'r ni {rW n2 {r')V{r - r')^ 3 {r')^ ni {r). (17) 

Let us perform the Bogoliubov transformation: 

b v = ^(iweini + v nv a' nl + u^ n a n 2 + v nv al 2 ), ( 18 ) 

which satisfies e„6„ = [b u ,H], leading to the corre- 
sponding Bogoliubov-de Gennes (BdG) equations. Such 
equations can be simplified by taking the more appro- 
priate variables A va = u„„ + u VTl , B VYl = v vn + v UYl , 
C vn = w„ n - u vri , D„ n = v VYl - v un . Then, the BdG 
equations take the form: 

e v A un = hcj n A vn , (19a) 
( v B vxl = -hw n B vn , (19b) 

+2N 1 d 2 {3(n', 0; n, Q)C vn , ~ g(0, 0; n, n')D vn ,}, (19c) 

+2N 1 d 2 {fl(n, n'; 0, 0)C„ n - g(n, 0; n', 0)IW}. (19d) 

In the following, we omit for simplicity v in the nota- 
tion for the A, B, C and D coefficients. We observe 
that, like in the homogeneous case, we have two dif- 
ferent branches of quasi-particles, one of them being a 
soft mode. The last two equations are the BdG equa- 
tions for the non-soft mode. From the properties of 
spherical harmonics and the 3j-symbols, it is possible to 
obtain that g(0,0;nlm,n'l'm') — g(nlm,n' I'm'; 0,0) = 
g (0,0; nlm, n'l' — m)5 m i- m g(nlm,0;n'l'm' ,0) = 
g{n'l'm', 0; nlm, 0) = g(0, 0; nlm, n'l' — m)8 m i tm . Also, 
V = 1,1 ±2 must be fulfilled. Employing these prop- 
erties, and changing to the variables F n \ m = C n i m — 
(-l)' Tl D„i_ m , Gnim = C n i m + {-T) m D n i-m-, one obtains 
that G„im = j-^F n i m , and 

^AN^hujr, 5(0, 0; nlm, n'l' - m)F n , Vm . (20) 

One can write: 

g(0,0;nlm,n'l'm') — 

{ ~ 1)(l '~ l)/2 2^ h(n ' l '' l ' ){l ' m]l ' m|20> ' (21) 



where do = y/K/2mu). In Eq. ( pi| ) the angular contribu- 
tion is given by 



(l',m\lm\20) 

I' 1 2 \ 



(-iy 



-mmO 



I' 12 
000 



5(2J + 1)(2Z'+1) 



nl/2 



47T 



(22) 



whereas the coefficients h(n, I; n', I') constitute the radial 
contribution and are of the form: 



h{n, I] n , I ) 



1 



3 ^5 2(«+«')/2 

r ( " + ^'+ 3 ) 



(^±3) r r (=^±a) r (^) ' 



(23) 



Then the problem reduces to finding the eigenvalues of 
the following equation: 

e 2 -Fnim = n 2 F n i m + \f2C,n 



V=l,l±2 n'>l 



J2(-l)^- l '^ 2 h{nl;n'l')(l'm\lm\20)F nn , m , (24) 



where ( is defined in Eq. ([To|), and e is in units of (?kc>) 2 . 
From the analysis of the spectrum one can observe that 
for £ > ( cr — 6.7 the energy of the lowest state of the 
subspace n = 2, m = becomes imaginary. Therefore, 
although the ground-state of the mixture of antiparal- 
lcl dipolar components is that of an ideal gas, the sys- 
tem will eventually become unstable for 7/ > rj cr due 
to the appearance of imaginary excitation frequencies. 
The value of £ cr compares very well with the numerically 
found value of Ccr — 7.3 for the onset of the instability. 

Finally, it is interesting to analyze the situation in 
which ( C 1, since in this situation the states with dif- 
ferent n do not mix, and one can analytically diagonalize 
in each n subspace. The first excited states are (in units 
of (huj) 2 ): 

e 2 (n = 1, 1 = 1, m = 0) = 1 + 4V27^C/15, 

e 2 (n = 1,1 = l,m = ±1) = 1 - 2 V /27^C/15, 

e 2 (n = 2,1= {2,0},m = 0) = 4 + (2 ± a/102) a/^C/^I, 

e 2 {n = 2, / = 2, m = ±1) = 4 + 2 ^2/^(/21. 

From these expressions it becomes clear that, for exam- 
ple, in the subspace of states with n = 2, m = (which 
is an admixture of the noninteracting quadrupole and 
monopole modes) even for £ ~ 0.1 deviations of more 
than 1% are expected. One can also observe that vor- 
tices with vorticity in the dipole direction are less ener- 
getic than those with vorticity in the plane perpendicular 
to the dipole direction. 



VII. CONCLUSIONS 

In this paper we have analyzed the ground state and 
the excitation spectrum of single- and two-component 
dipolar condensates. For the case of single-component 



10 



BECs we have shown that their stability properties are 
determined by the trapping geometry. In particular, 
for sufficiently pancake-like traps the condensate is al- 
ways stable independent of the number of particles. We 
have then analyzed the excitation spectrum of the single- 
component condensate by analyzing the response of the 
condensate to small perturbations and comparing the re- 
sults with analytical calculations based on a variational 
approach. We have discussed in detail the nature of the 
instability and associated it with vanishing of frequency 
of one of the excitation modes. The scaling behavior 
of this frequency was also analyzed. We have discussed 
different possible scenarios in which the dipole-dipole ef- 
fects could have observable effects even in the presence of 
a dominant short-range interaction. This analysis could 
be of special interest for ongoing experiments on atoms 
with large magnetic moments, such as chromium [[39). 
We have provided guidelines for experimental parame- 
ters corresponding to largest discrepancies between the 
cases with and without dipole-dipole interactions. In the 
second part of the paper we have analyzed the properties 
of a two-component BEC of dipolar particles. In particu- 
lar, we have studied the stability of the ground state as a 
function of the relative density of both species as well as 
the appearance of phase separation in the binary conden- 
sate. Finally, we have obtained the excitation spectrum 
for this particular physical system and discussed the na- 
ture of its instability. 
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